2.11 ATOMIC HYDROGEN DISTRIBUTION 


N. Tabarie 


Introduction 


Molecular hydrogen (Ha) has been identified in Titan's atmosphere by its 
quadrupole absorption line, with an estimated total quantity of 5 km-A (Trafton 
1972). There must then be some atomic hydrogen (H) in Titan's atmosphere, which 
may be accessible to observation in Lyman-a at 1216 A. 

In this study we have examined several possible vertical distributions 
with the constraint of 5 km-A as a total quantity, and have calculated, with 
approximations, the corresponding vertical distributions of atomic hydrogen. 

It was found that the H distribution is quite sensitive to two other parameters 
of Titan's atmosphere: the temperature and the presence of other constituents. 

The escape fluxes of H and H 2 were also estimated as well as the consequent 
distributions trapped in the Saturnian system. 


Models of H 2 Distribution 


The constraint of 5 km-A of H 2 makes impossible an atmosphere of pure H 2 , 
since it would escape so fast that the total mass of Titan would vanish in less 
than 10® years. Additional constituents with molecular weight higher than H 2 
have to be present to accommodate the 5 km-A of H 2 . Two constituents were con- 
sidered separately; (1) N 2 after a suggestion by Hunten (1972),;with various 
mixing ratios ranging from 10”^ to 10“^ for H 2 /N 2 ; and (2) CH 4 , which was posi- 
tively identified by Kuiper (1944), with a mixing ratio of 1. 

For the sake of simplicity the atmosphere was considered to be spherically 
symmetrical and isothermal, with a temperature of either 80, 100, or 120°K. By 
analogy with the terrestrial atmosphere, three different zones were considered: 
The turbosphere where the mixing is constant, the diffusosphere where each con- 
stituent follows its own scale height, and the exosphere where collisions are 
not important. The turbopause level was defined as the altitude where the 
total number density is lO'-^ mols cm"^. The exobase level was defined as the 
altitude Z^. where the density of H 2 is 10 ^ mols cm"^^ except if it was found to 
be higher than the blow off level (kinetic energy equal to gravitational energy) , 
in which case the exobase was located at the blow off level. These levels are 
indicated in Table 2-8, together with the escape fluxes of H 2 at the exobase 
estimated from the Jeans formula. The various calculated H 2 distributions are 
indicated in Figure 2-39. The main feature of these H 2 distributions is the 
existence of a very large diffusosphere when compared to the terrestrial case. 
This is due to the low gravity of Titan and to the high total quantity of H 2 . 


Models for H Vertical Distribution 

Photodissociation of H 2 and CH 4 by the solar flux was considered as the 
only source of H, compensated by escape at the exobase and by three-body re- 
combinations at lower altitudes. For each H 2 profile a corresponding H profile 
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Table 2-8a. H 2 Fluxes for Titan Atmospheric Models 


TEMP. 

MIXING 

RATIO 

TURBO- 

PAUSE 

Is. 

km 

EXOBASE 

Zc 

km 

DENSITY 
AT GROUND 
{H 2 ) 

mols cm- 3 

DENSITY 
AT Zc 
{H 2 ) 

mols cm 3 

EFFUSION 
VELOCITY 
(H 2 ) , 
cm sec 1 

ESCAPE FLUX 
(H 2 ) 

mols cm- 2 sec 1 

TOTAL FLUX 
(H 2 ) 

mols sec 1 

FLUX 

AT GROUND 
(H 2 ) 

mols cm- 2 sec-i 


10-3 

520 

1570 

11.0 X 1018 

102 

2.3 X 102 

2.3 X 103 

4.8 X 1022 

5.9 X 109 

80°K 

(H 2 -N 2 ) 

10-2 

470 

3420 

10.5 X 1018 

102 

1.3 X 103 

1.3 X 1018 

5.9 X 1028 

7.0 X 1018 


io-‘ . 

460 

8950 

9.5 X 1018 

10 ? 

7.0 X 103 

7.0 X 101“ 

1.6 X 1038 

2.0 X 1012 


10-3 

670 

2410 

8.0 X 1Q18 

102 

1.5 X 103 

1.5 X 1018 

4.6 X 1028 

5.6 X 101“ 

100“K 

{H 2 -N 2 ) 

10-2 

610 

6050 

8.0 X 1018 

102 

6.8 X 103 

• 6.8 X 1018 

6.3 X 1028 

7.7 X loll 


10-1 

590 

12060* 

7.0 X 1018 

3.6 X 102 

1.4 X lO" 

5.0 X loll 

1.3 X 1031 

1.6 X 1Q13 


10-3 

840 

3470 

6.4 X 1018 

■ 102 

5.5 X 103 

5.5 X 1018 

2,4 X 1028 

3.0 X 1011 

120°K 

(H 2 -N 2 ) 

10-2 

750 

9600* . 

6.3 X 1018 

1.7 X 102 

1.6 X 10" 

2.6 X 1011 

4.8 X 1038 

5.8 X 1012 


10-1 

730 

9600* 

5.6 X 1018 

1.5 X 108 

1.6 X 10 " 

2.3 X 1012 

4.2 X 1031 

5.1 X 1013 

SO^K 

(H 2 -CH 4 ) 

■ 

1470 

15710* 

2.7 X 1018 

2.7 X 108 

1.2 X 10“ 

3.2 X 1012 

1.4 X 1032 

1.7 X 101“* 

100 °K 

{H 2 -CH 4 ) 

■ 

2070 

12060* 

2.1 X 1018 

2.3 X 108 

1.6 X lO" 

3.2 X 1013 

8.6 X 1032 

101 ? 


* Exobase * Blow-off level in this case. 
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Table, 2-8b„, ,,.H Fluxes for Titan, Atmospheric Models,. 


^ . 

. 



. . 




TEMP. , 

'mixing . 

.RATIO 

TURBO- 
PAUSE 
. Za 
kjl 

EXOBASE 

Zc 

km 

BENSITY 

AlZc 

(H) 

atoms cm" 3. 

ESCAPE .flux ' , 
atoms cm ‘sec ' , 

TOTAL FLUX 
■(H) ■ 

atoms sec * ' 

FLUX RATIO 
H/H 2 



52d 

, 1570 

j 7.50 X 10“ 

3.7, X 103 

8.0, X 1026 

1.6 X 10-3 

80°K 

- (H 2 -N 2 ) . 


470 

3420 

, 2.50 X 10“ 

2.7 X 10® 

1.2 ,x 1022 

2.0 X 10-2 



460 

8950 

7.80 X 103 

1.7 X 103 

2.8 X 1027 

1.7 X 10-3 


10-3 

. 670 

2410 

, 2.30 X 10“ 

2.9, X 10® 

9.0 X 1028 

2.0 X 10-2 

100°K 

(H 2 -N 2 ) 

10-2 , , 

610: 

, ,6050 - 

7.40 X 103 

1.7 x,108 

1,5 X 1027 

2.0 X 10-3 


10-1 

590 

12060* 

1 10“ 

3.0 X 108 

8.0 X 1,02 7 

■ 

6.0 X 10-“ 


10"3 

. 840 

3470 

1.10 X 10“ 

2.5 X. 10® 

1.1 X 1027 

6.5 X 10'3 

120°K 

(H 2 -N 2 ) 

10-2 

750 

9600* 

. . 5.60 X 103 

1.8 X 108 

3.3.x 1027 

■ 6.8 X 10'“ “ 
'1 

2.0 X 10-“ 


10-1 

730 

- 9600* 

1.60 X 10“ 

5.-5 X 10« 

1.0 X 1028 

SO^K 

(H 2 -CH 4 ) 

1, 

1470 

15710* 

3.00 X 10“ 

0 

X 

0 

CO 

3.3 X 1028 

2.0 X 10-“ 

100°K , 
(H 2 -CH 4 ) 

ll 

2070 

12060* 

' 1.25 X 103 

3.7 X 10^ 



1027 

10-8 


* Exobase = Blow-off level in this case. 
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Figure 2-39. Profile of molecular hydrogen for H 2 -N 2 and H 2 -CH 4 model atmos- 
pheres at temperatures of 80°K and 100°K. is the turbopause 
level and Z^, is the exobasei level. 
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was calculated along the vertical at the sub-solar point by solving simultane- 
ously the diffusion equation and the continuity equation in the spherical case: 


= -D(r) 


dn 

dr 


H(r) 


d((ij^Cr) 

dr 


= P(r) 


L(r) - I (r) 


where: H(r) 

Hg(r) 

PW 

h(r) 

D(r) 


scale height of atomic hydrogen, 
scale height of the atmosphere. 

the diffusion flux of atomic hydrogen in ats cm"^ sec"^, 
production rate of H through photodissociation of H 2 , 
loss rate of H by three-body recombinations with H 2 and N 2 , 
diffusive coefficient: 

2 X 10^^ 

D(r) = in the case of diffusion of H in an 

"H2(r) 

atmosphere of H 2 , 

3.19 X IQlS H Cr) 

D'(r) = : ^ in the case of diffusion of H, 

”N2(r) 

in an atmosphere essentially of N 2 . 


The exospheric distribution wqs calculated with no satellite particles from 
Chamberlain's theory. 

For H 2 -N 2 model atmospheres the effect of the mixing ratio CH 2 /N 2 ) at 
ground level is illustrated in Figure 2-40 for T = 80°K. When this ratio in- 
creases the density maximum Of H, located in the turbosphere, decreases, while 
the density of H in the diffusosphere increases substantially. Above lO**^ km 
altitude in the exosphere, the density is insensitive to the mixing ratio. 

The effect of the isothermal temperature (T) on the H density is illus- 
trated for an H 2 ~N 2 mixing ratio of 10"^ in Figure 2-41. The H density increases 
substantially with T in the diffusosphere. Again the density in the exosphere is 
not very sensitive to the temperature. 
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10 ^ 10 ^ 10 ^ 10 ® 


ALTITUDE, Z (KM) 

Figure 2-41. Profile of atomic hydrogen for H 2 -N 2 model atmospheres for 

temperatures of 80“K, 100“K, and 120°K. A mixing ratio, H 2 /N 2 
of 10" 1 IS assumed. is the turbopause level and is the 
exobase level. 



For H2-CH4 model atmospheres the effect of temperature is illustrated in 
Figure 2-42 for a mixing ratio of unity. The density profiles of H are quite 
different than for. the H2-N2 models; the density at 5 x 10^ km and on into the 
exosphere is much higher for the same temperature of 80°K. The effect of T is 
at variance with H2-N2 models in that an increase in T results in a decrease 
in the H density in an H2-CH4 model. This difference between the two types of 
models is due partly to the higher mixing ratio in the H2-CH4 model and partly 
to the increase of H production through photodissociation of CH4, This later 
process is imderestimated in our approach, since only photodissociation giving 
CH3 + H was considered. 


H and H2 Escape Fluxes 


The total fluxes of H and H2 indicated in Table 2-8 were calculated by 
assuming spherical symmetry which might lead to an overestimation by a factor 
of 2 for H. In any case, the total flux of H is found to vary quite widely 
with the atmospheric parameters between 7.9 x 10^® atoms sec”^ for the H2-N2 
model to 3.3 x 10^® atoms sec"^ for the H2-CH4 model at 80°K. The flux of H2 
is even larger than the flux of H for all models by a factor of between 6 and 
10®. For H2-CH4 models the H2 flux is much larger than for H2-N2 models, and 
a hydrodynamic approach would be more appropriate than the assumption of hydro- 
static equilibrium. 

It has been suggested by Brice and McDonough (1973) that a large escape 
flux of H2 from Titan would result in a toroid of hydrogen around Saturn. 
Dennefeld (private communication) has calculated the distribution of atomic 
H resulting from the escape of H and H2. He found that a flux of 2 x 10^® 

atoms cm"® sec"^ of H would result in a mean density of 22 ats cm"® in the 

whole volume indicated in Figure 2-43 around the orbit of Titan. With a den- 
sity of 22 ats cm"®, such a toroid of H should be easily detected by its reso- 
nant Lyman-a emission, with an intensity of several hundred Rayleighs. In 
addition, since 10% of escaping H2 molecules will be photodissociated into 2 H 

atoms when they are around Saturn, a flux of 10®l H2 molecules (corresponding 

roughly to all encountered blow off conditions) would result in a much larger 
mean density of 22 x 100 = 2200 ats cm"®, greatly increasing detectability of 
the toroid. 


Summary 


The distributions which we have obtained show that a determination of the 
vertical distribution of atomic hydrogen could easily be fitted to a model and 
that the two significant parameters, temperature and mixing ratio, could then 
be determined. However, since the differences between H2-N2 and H2-CH4 models 
are not very large, a high spatial resolution in the measurements would be neces- 
sary. 


Hun ten : Going back to my discussions of the limiting escape flux, the distribu- 

tions of H and H2 are independent of K (the eddy diffusion coefficient) until it 
reaches a value even larger than was assumed here (2 x lO® cm® sec"^). For very 
large values, the densities are decreased. (K is related to the number density 
at the turbopause, Ug^, by K = 2 x 10l®/ng.) What is the storage time in the 
toroid for hydrogen atoms and molecules? 
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Figure 2-43. Cross-section o£ toroidal zone H and H 2 concentrations in the 
Saturnian system. 
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Tabarle : About 2 x 10® sec which is determined by charge exchange and photo- 

ionization of the atoms. 


McDonough : You will be able to see Lyman-a both from resonance scattering by 

H and photodissociation of H 2 . Do you have any idea what the ratio is? 


Hun ten : This was an issue after the Mariner 5 flyby of Venus. Equal intensities 
from the two sources required an H 2 /H ratio of about 10 ®; 


Blamont : ...which may be what we will have on Titan and in the toroid, but the 
line from H 2 fluorescence will be much wider and a good measurement of the line 
shape will allow you to discriminate. 


Rasool : What will the Mariner Jupiter/Satum (MJS) mission be able to do near 

Titan? 


Broadfoot : MJS will have two spectrometers, one observing airglow, with a 6 

arc-minute field, and one observing solar occultation (if it goes through Titan's 
shadow) with a 1-minute field. Occultation will give H, H 2 , and CH 4 . • Airglow 
will give H, but with poorer vertical resolution. 


Blamont : It seems doubtful that MJS can detect H by occultation, because the 

solar Lyman-a line is so much wider than the planetary absorption. 


Hun ten : The height resolution is not too bad if you can get close enough to 

Titan. One arc-minute at 30,000 km is 10 km, which isn't bad at all. 


Sagan : If your instrument can measure as faint as 100 Rayleighs, it, would be 

exciting to map the glow from the toroid or tail or whatever is there. 


Broadfoot : We will be able to do that as we approach the planet; we can take 

days to scan over the whole Saturnian system. 



